Effect of gap suppression by superfluid current on nonlinear microwave response of 

d-wave superconductors 
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Recently several works have focused on the intrinsic nonlinear current in passive microwave filters 
as a tool for identifying d-wave order parameter symmetry in the high T c cuprates. Evidence has 
been found for d-wave pairing in YBCO and further work has ensued. Most of the theoretical work 
has been limited to low temperature because it has not included the effect of the superfluid current 
on the energy gap. We find that this effect leads to important corrections above T ~ 0.2T C , while 
leaving the 1/T low temperature behavior intact. A twofold increase in the nonlinear coefficient at 
temperatures of order ~ 0.75T C is found, and as T — > T c the nonlinearity comes entirely from the 
effects of the superfluid current on the gap. Impurity scattering has been included and, in addition, 
signatures for the case of d+s-wave are presented. 

PACS numbers: 74.20.Rp,74.25.Nf,74.72.Bk 



I. INTRODUCTION 

Interest in generating novel methods for probing the 
order parameter symmetry in superconductors has driven 
the development of new techniques and experimental con- 
figurations, such as angle-resolved electron tunneling^] 
or angle-resolved magnetic field dependence of specific 
heatQ, which can provide unambiguous signatures of or- 
der parameter nodes and their location. One such pro- 
posal is associated with the measurement of the intrinsic 
nonlinear current in passive microwave filters. Such de- 
vices give rise to to third order intermodulation effects 
which, while detrimental for practical applications in 
superconducting communication filter technology, make 
these devices ideal for examining issues of order parame- 
ter symmetry. Excellent progress in the field of high tem- 
perature superconductivity has been made in this area. 

Initially, Yip and Sauls 0,3 proposed the examination 
of the nonlinear Meissner effect (NLME) for evidence of 
d-wave gap symmetry in the cuprates. Their predictions 
were not confirmed by experiments at the time 3, how- 
ever, Dahm and Scalapino (DS)@,0| proposed to exam- 
ine a related quantity: the intermodulation current or 
distortion (IMD) which arises from the nonlinear induc- 
tance resulting from a quadratic dependence of the pen- 
etration depth on the superfluid current. In their work, 
a signature of the d-wave gap would be found in an up- 
turn in the temperature dependence of the nonlinear co- 
efficient at low temperatures, in contrast to exponential 
decay for an s-wave gap. Such evidence of an upturn has 
been found in YBCO films by several groups]!, 03- E3 
with excellent agreement with the DS theory. This suc- 
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cess has led to further calculations to examine the issue 
of nonlocal effects for this quantity^]] (thought to be 
the possible reason for not seeing the NLME proposed 
by Yip and Sauls ^3) an< ^ the d-wave signature remains 
robust. More recently, measurements of the nonlinear 
current have been made near T c on a range of different 
films and the data has been analyzed in terms of the DS 
approach with impurity scattering being used to explain 
variations between the films|l3|. However, in this work 
and others^l ^| which extrapolate the DS theory to 
high temperature, the effect of the superfluid current on 
the gap has been neglected as an approximation. We 
show here that inclusion of this effect has significant im- 
pact on the nonlinear coefficient for T > 0.2T C , requir- 
ing a reanalysis of previous results. Indeed, when the 
current-dependence of the gap is considered, the nonlin- 
ear coefficient shows a twofold increase at T ~ 0.75T C 
over the value calculated using the approximation of a 
current-independent gap. 

Finally, we note that other theoretical works have ex- 
amined the intrinsic nonlinear current for two-band su- 
perconductors and M g B9 |15lll6j | and for one-band s-wave 
superconductors P. Hq|. Effects of the superfluid on the 
gap have been included in Refs. P.ll6|. 

Our paper is structured as follows: in the next section, 
we summarize our theoretical approach which allows for 
both the inclusion of the current dependence in the gap 
and impurity scattering (from unitary to Born limit). 
Strong electron-boson coupling effects are also available 
in this formalism. In Section III, we discuss the results il- 
lustrating the corrections to the simplified theory at finite 
temperature and we revisit the experimental situation, 
including the issue of impurities. We end by providing 
predictions for an admixture of d- and s-wave symmetry 
as has been recently suggested by angle-resolved electron 
tunneling experiments on YBCO0. We form our con- 
clusions briefly in Section IV. 
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II. THEORY 



In this work, we evaluate the full current from the stan- 
dard expression given for the imaginary axis Matsubara 
representation and modified for a d-wave order parameter 
A(0) = Acos(2(9) in two-dimensions: 0, H [H GjJ 
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This expression contains both the condensate current and 
the quasiparticle current due to excitations.^ The an- 
gle a measures the direction of the current with respect 
to the order parameter antinode, with a — indicat- 
ing the current in the antinodal direction and a = 7r/4 
for the nodal direction. The other notation is standard 
with e the electric charge, m the electron mass, T the 
temperature, n the electron density and vf the Fermi 
velocity. The superfluid momentum q s enters through 
s = Vpq s , and, importantly for the results of this paper, 
it also enters the equations for the Matsubara gaps and 
renormalized frequencies and, as a result, the current will 
decay the gap. The equations for the Matsubara gaps 
A„ = Z n A n and renormalized frequencies Q n = Lu n Z n 
modified for d-wave symmetry are: 
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where the Matsubara frequencies are u) n — irT(2n — 1), 
for integer n, and we have included the possibility of 
impurity scattering from the unitary (c = 0) to Born 
(c — » oo) limit via the last term in Eq. J2J| which requires 
self-consistency through il n . Here, T + is proportional to 
the impurity scattering rate and c is related to the the 
s-wave scattering phase shift. 0] Note that an impurity 
term does not appear in Eq. J5J in d-wave as it does in 
s-wave. This is because in d-wave it averages to zero. 
Finally, in general, the kernel of these equations would 
normally be based on a momentum-dependent electron- 
boson spectrum. To mimic this unknown spectrum in 



an approximate manner, a parameter g is introduced 
to represent that the interaction in the w-channel could 
be different from that in the A-channel in the case of a 
momentum-dependent interaction that would give rise to 
a d-wave order parameter symmetry. Likewise, there is 
no cos(20) factor in the numerator of the w-channel as 
there is in the A-channel reflecting that the interaction 
in the renormalization channel is taken to be isotropic 
to first order. pH |22T| Thus, the electron-boson spectral 
function, which we denote by a 2 F(fl), enters X(n — m) 
as follows: 
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Here, to take the limit of these equations to give the stan- 
dard BCS result for d-wave in our numerical evaluation, 
we take the electron-boson spectrum to be a delta func- 
tion at high frequency. Likewise, we exclude renormaliza- 
tion effects which are not based on impurities by taking 
g to be zero. This gives the BCS gap ratio in d-wave to 
be 2Ao/kT c = 4.28. When we wish to consider strong- 
coupling effects corresponding to 2A / kT c = 5, for ex- 
ample, then we take g to be a finite value and use a delta 
function at lower frequency for the boson spectrum p2|. 

To extract the nonlinear coefficient that is relevant to 
the passive microwave filters and hence can be used as 
a sensitive probe of order parameter symmetry, we can 
assume that j s can be expanded to third order for small 
q s as: 
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=fl(^)-«r,( 



q s VF 

An 



3n 



(6) 



where jo — neA^/ (mvp). Note for strong coupling, we 
replace we replace the q s VF in this formula by q s v f/ Q + 
A), where A is the mass renormalization parameter 
From this Dahm and Scalapino define 



b(T) 



P(T) 



[n 8 (T)/n] 3 ' 



(7) 



the square of which is related to the third order intermod- 
ulation power in microwave filters. Hence, measuring the 
intermodulation power provides a measure of the nonlin- 
ear coefficient. As we calculate the full q s dependence of 
j s , using Eqs. Q-© with no approximation of taking 
the gap to be independent of q s , as was done in previ- 
ous works, it is easiest to extract these quantities directly 
from our numerical data. To do this we form the quan- 
tity of js/ls versus q 2 which is a straight line at low q s 
and from this we obtain the superfluid density from the 
intercept and the nonlinear coefficient from the slope. 
By this method, we have confirmed previous results at 
low temperature and can proceed to examine the issue of 
higher temperatures where the current reduces the gap 
even at low q s . We have also demonstrated this method 
for one-band and two-band s-wave superconductors. |l6j 
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III. RESULTS 

In Fig. ^ we show the current as a function of q s in 
the two major directions, along the node and antinode, 
at both low and high temperatures. This was done using 
the equations above and illustrates that we can reproduce 
correctly the T = results in the literature jU El and 
that we can also evaluate the current at high tempera- 
ture in this formalism. The Matsubara formalism is also 
ideally suited for including impurity scattering. There 
are few points to note in this figure. The q s in the order 
parameter is essential to obtain these curves and it is the 
decay of the order parameter by the superfluid current 
that causes the current j s to drop dramatically beyond 
the peak (otherwise, if A(q s = 0) is used, the curves de- 
cay slowly to zero as q s — > oo, for example, as l/g s for 
T = BCS s-wave). Furthermore, the effect of the q s in 
the gap becomes even more important at low q s when the 
temperature is approaching T c . Finally, for high temper- 
atures near T c , the current is fairly independent of the 
angle a. 




0.0 0.2 0.4 0.6 0.8 1.0 



q s v F/ A o 

FIG. 1: The normalized current j s /jo as a function of 
QsWf/Ao, where jo = neAo/(mvF) and Ao is the energy gap 
at T — 0. Shown are the low temperature BCS curves for two 
directions: the current in the antinodal direction with a = 
(dashed line) and the nodal direction with a = 7r/4 (solid), 
given for a reduced temperature t = T/T c — 0.1. The in- 
set shows that the curves overlap for T near T c (in this case, 
t = 0.95). 

This latter feature is seen more clearly in Fig. [21 where 
we show the nonlinear coefficient /3(T) for d-wave in the 
clean limit as a function of temperature for the two di- 
rections just discussed (solid and short-dashed curves). 
Once again, it is seen that while the two curves are dif- 
ferent at low T, as T — > T c , the anisotropy is reduced 
and disappears at T c . Indeed, we can obtain an analyt- 
ical value for (3(T) at T = T c which is 0.651 and this is 
confirmed by the numerics shown in Fig. EI It is a ls° i m ~ 



portant to note that the same procedure which gives rise 
to these curves also provides the normalized superfluid 
density which is shown as the solid curve in the inset 
of Fig. |21 This superfluid density curve is independent 
of the direction of the current as expected as it arises 
from the q s — > limit and this curve is exactly the same 
as that which is calculated by standard formulas for the 
penetration depth. Here it is extracted from our j s ver- 
sus q s curve as explained in the theory section, verifying 
the accuracy of our method. 
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FIG. 2: The nonlinear coefficient /3(T) as a function of T/T c , 
shown for two directions: a — n/4 (solid line) and a = 
(short-dashed). The dotted curves are for the same two di- 
rections but with the approximation of neglecting the q a de- 
pendence in the gap, i.e. A(q a = 0). The inset illustrates 
via a log-log plot that the low temperature behavior varies as 
Ao/24T for a = and Ao/12T for a = ty/4 (these expressions 
are shown as the long-dashed lines in both cases). 

Also shown in Fig.|21by dotted line type are the curves 
for P(T) for a = and 7r/4 when the approximation of 
A.(q s = 0) is taken (i.e., the s is set equal to zero in 
Eqs. so that the gap is not modified by the super- 

fluid current). This is a central point of our paper, that 
while this approximation works well at low temperatures, 
one sees from the comparison of the dotted curves with 
their respective short-dashed and solid ones, that this 
approximation breaks down for T > 0.2T C and produces 
significant deviations as T — > T c . Indeed, from physical 
grounds one does not expect the nonlinear current to go 
to zero as T — > T c , but rather to increase. From the point 
of view of device applications which would typically op- 
erate at about 0.5T C or higher, this effect of decaying the 
gap by the superfluid can introduce a factor of 1.5-2 in- 
crease in the nonlinear response of the device. Finally, we 
confirm in Fig. |21 that the original use of this approxima- 
tion for low temperatures is robust and the 1 /T signature 
of the d-wave gap discussed by Dahm and Scalapino0,0 
remains intact. We find by our procedure, the same result 
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as determined by Dahm and Scalapino analytically, that 
P(T,a = 0) ~ A /24T and (3(T,a = vr/4) ~ A /12T 
for T — > 0, which is illustrated by the log- log plot in the 
insert of Fig. [21 




o.o I I 
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FIG. 3: The nonlinear coefficient 0(T) versus T/T c q for vary- 
ing impurity scattering in the unitary limit (c = 0) and with 
a = 7r/4. The inset shows the corresponding curves for the 
superfluid density n s (T)/n. Curves are given for: the pure 
limit (solid), F + /T c0 = 0.0082 (short-dashed), 0.0404 (long- 
dashed), 0.0697 (dot-short-dashed), 0.1432 (dot-long-dashed). 

Turning to the case of impurities which were discussed 
previously by Dahm and ScalapinoJfiJ and Andersen et 
aL|l3|. we comment on the modifications that occur due 
to the current in the gap. Again, the results of Dahm 
and Scalapino which are given for T < 0.2T C remain ro- 
bust, however, the results of Andersen et al. which focus 
on high temperature are necessarily modified when the 
approximation of a (^-independent order parameter is re- 
moved. This is shown in Fig. [21 where we show both (3(T) 
and n s (T)/n for varying impurity content (here we use 
the unitary limit for illustration, as was done in Ref. [l3|). 
The (3(T) curves are shown only for a = n/4. Once again, 
for T > 0.2T C , deviations occur in the manner discussed 
before. The significant feature is that the impurity scat- 
tering reduces the T c relative to the pure case which has 
T c o and hence for fixed temperature of T ~ 0.75T c o, for 
example, with varying impurity content, the decay of the 
gap by the superfluid current has an even greater effect 
for increasing impurity content and so it must not be 
neglected. 

We now turn to a re-examination of the comparison of 
theory with some of the data that exists in the literature, 
maintaining the same analysis that was done previously. 
First, we consider the clean limit YBCO data presented 
by Oates and coworkers^. In Fig 01 we reproduce this 
data for the normalized IMD along with our calculations 
for BCS d-wave, where 2A = 4.28fcT c . The IMD power 
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FIG. 4: Plot of \og(Ab 2 av {T)) versus log(T/T c ) for s-wave 
(dashed curve), BCS d-wave with 2A /fcT c = 4.28 (solid), 
and strong-coupling (SC) d-wave with 2Ao/kT c = 5 (dotted). 
The curves are compared with the experimental data (solid 
dots) taken from Ref. [9j. The inset shows the same curves 
for b(T) versus T/T c . 



has been shown to be related to b 2 (T), i.e. 

P IMD cx b 2 (T) (8) 

and hence we plot b 2 (T) in the figure and use an ad- 
justable factor A to match the data with the theory as 
was done by Oates et al.. Likewise, we have averaged 
the two directions for a — and a — 7r/4 after the man- 
ner of Oates et al. as their measurement averages over 
all directions. The low temperature part of this log- log 
plot does not change from previous d-wave calculations 
but the high temperature part is increased and the agree- 
ment with the data is not as good in this regime (we have 
chosen to fit the low temperature part of the curve to the 
data). However, compared to the s-wave calculation, the 
evidence in support of d-wave remains striking. Oates 
et al. obtained a better fit with the d-wave theory be- 
cause Dahm and Scalapino used a larger gap ratio in their 
BCS theory which is supported by other experiments. To 
illustrate this within our formalism, we can do a strong- 
coupling calculation where we take g = 1 in Eq. [21 and 
move the electron-boson spectral function to lower fre- 
quency. In this manner, we can obtain 2A /A:T C — 5 in 
d-wave corresponding to a mass renormalization param- 
eter A = 5. In this case, the result (shown as the dotted 
curve and plotted with a new A to fit the low T region) 
does indeed move toward a better fit with the data over 
the full temperature range. As a technical note, to obtain 
a gap ratio of 2Ao = 6fcT c in this Eliashberg formalism, 
we would have to take A to extremely large values and we 
would then be in the asymptotic regime which limits the 
value of the gap ratio to go no higher than about 6.5. [22^ 
A much better fit is obtained with a gap ratio of about 6 
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but in our model the A is unphysical. To overcome this, 
in order to produce higher gap ratios with more physical 
values of A, would require a spin-fluctuation calculation of 
the type given in Ref. pEj , which would include feedback 
on the electron-boson spectral density itself, which sup- 
presses it at low lo and introduces coupling to a new reso- 
nant mode that grows in amplitude as T is reduced. This 
is beyond the scope of this work. Other work which has 
discussed this data has been that of Agassi and Oates[l"H 
where they considered nonlocal effects, however, the cal- 
culation is a BCS one with no inclusion of the effect of 
the superfluid current on the gap and so it is difficult 
to say how both this and strong coupling might modify 
their results. 

Finally in Fig. in the inset, we give for reference 
the unsealed b(T) versus T curves. The 1/T divergence 
is clear in comparison with the exponential decay of the 
s-wave case and the strong-coupling calculation is above 
the BCS one as is expected from previous work on s-wave 
superconductors . 
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FIG. 5: Plot oilos(j 2 (T)/j c ) versus A(T)/A for T = 0.75T c0 . 
Andersen et aZ.|l3f have defined a quantity (jb/ic) 2 which is 
equal to our l/b(T). The data (open circles) is taken from 
Ref. The solid curve is for BCS d=wave with A(q s ) and 
the dashed is for the case where the approximation A(q s = 0) 
has been used. The inset shows the reduction in T c from the 
pure value of T c o that is implied by the values of A(T) at 
T = 0.75T c o normalized to the pure value for T = 0, Ao. 

We revisit the issue of impurities in Fig. El This fig- 
ure is based on a similar one presented by Andersen et 
aL^j|. Several films were examined at 75K and the non- 
linear coefficient was found to vary from sample to sam- 
ple, as did also the penetration depth. The analysis in 
the paper used the DS theory with impurity scattering 
to argue that the data could be understood by assuming 
varying impurity content from one film to another. We 
re-examine this issue because at this high temperature 
the pairbreaking effect of the superfluid current on the 



gap should cause significant chan ges (our Fig. [3] should 
be contrasted with Fig. 3 of Ref. |R|) an( i this was n °t 
included in the original analysis. In Fig. [S] we show our 
calculation at T = 0.75T c o for the variation of the non- 
linear coefficient with respect to the variation of the pen- 
etration depth when unitary scattering is included (like 
Andersen et at, we have checked Born scattering and find 
there is essentially no difference in the result presented 
here). We have also taken a — ir/4 as at this temper- 
ature the anisotropy is almost completely gone, which 
we have checked. It should be noted that, once again, 
there is an adjustable parameter j c which can be used to 
scale the theory to overlap with the data. Andersen et 
al. kept the theory fixed and used a value of j c to ad- 
just the data. Here we choose to keep the scaled data as 
it was presented in the original paper and so we adjust 
our theory by a scale factor to overlap with the data (if 
we had adjusted the data instead, this would have in- 
creased the j c value used for scaling the data by a factor 
of about 1.35). Doing so, we find with the solid curve 
that there is reasonable agreement between theory and 
data. However, we note, referring the reader to the inset 
of the figure, that the range of variation in penetration 
depth, if explained via impurity scattering, would i mply 
a reduction in T c of about 15-20% from the pure case|26| 
and yet the experimental data indicates a T c variation 
of only about 2%. This discrepancy was not noted in 
Ref. [13j and it does pose a problem for this interpreta- 
tion. However, continuing with this explanation, we have 
compared our result with the case of taking A(q s = 0), 
shown as the dashed curve and using the same j c scale 
factor, and we find that while some deviation may be cor- 
rected by choosing a different j c , which would shift the 
curve, the slope is also different. While the interpreta- 
tion of the data via impurity scattering may be open to 
question due to the large variation in T c required, more 
experimental data may help to reduce the scatter and 
provide further tests of the theory. 

Recently, a novel probe of angle-resolved electron tun- 
neling was used to determine the order parameter sym- 
metry in YBCO and this work has suggested that the 
order parameter is not pure d-wave but an admixture of 
d+s with about 15% of s-wave. This is not the first time 
that it has been suggested that YBCO may have a small 
s-wave component and it is important to assess the con- 
sequences for the intrinsic nonlinear current and whether 
there would be any signatures unique to a d+s-wave or- 
der parameter. As a result, we now briefly consider such 
a symmetry, where we assume that the s-wave compo- 
nent is small compared with the d-wave part, such that 
nodes still exist but they are shifted from 7r/4. Taking 
A(#) = A [a cos(26*)+c] , we evaluate analytically the main 
features expected in the limits of T — > T c and T — > 0. As 
T — > T c , we find that we can define a quantity related to 
b(T): 
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FIG. 6: The effect of a d+s-wave gap on the anisotropy near 
T c as measured via b(T c ) (upper frame) and at low tem- 
perature as measured by (3(T) normalized to Ao/24T (lower 
frame) . The anisotropy is plotted as a function of the current 
direction relative to the antinode, a. Curves are shown for s- 
wave (dotted), d-wave (short-dashed) and varying percentage 
of s- wave: 10% (dot-dashed), 15% (solid), 20% (long-dashed). 



which gives {T c / A ) 2 b(T c ) = 0.0599 for d-wave and and 
0.0266 for two-dimensional s-wave. In both cases, this 
quantity is isotropic with respect to the direction of the 
current, however, for d+s-wave, we have 
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and there arises an anisotropy as a function of a in the 
nonlinear coefficient as T — > T c , which is not there in pure 
d-wave. Indeed, it can be sizable as shown in Fig. (up- 
per frame) for different percentages of the s-wave compo- 
nent. For 15% s-wave, the anisotropy is 20%. Note also 
that the anisotropy is twofold. At low temperature, a d- 
wave order parameter already shows a fourfold anisotropy 
in (3(T) as a function of a. If we plot /3(T) normalized to 



Ao/24T as is shown in the lower frame of Fig- El we see 
that the fourfold d-wave anisotropy shifts to twofold for 
d+s-wave and also exhibits a shift in the position of the 
maxima (indicating a shift in the position of the nodes). 
The magnitude of the anisotropy, which was a factor of 
2 for d-wave, is now even greater in d+s-wave. The ana- 
lytic formula for this anisotropy, assuming c < a, is 



TO) - ^ 

24T 



1 — - cos(2a) 
a 



1- - 



sin 2 (2a) 



V~l {c/af 



(14) 



Thus, there should be observable signatures of an s-wave 
component should it be possible to do angle-resolved 
measurements of intrinsic nonlinear current. 



IV. CONCLUSIONS 

In summary, we have calculated the intrinsic nonlinear 
current of a d-wave superconductors, including the effect 
of the superfluid current on the order parameter. We 
find that the low temperature 1 /T behavior of the non- 
linear coefficient, proposed as a test for d-wave symmetry 
by Dahm and Scalapino, remains unchanged. However, 
at temperatures above about 0.2T C , the approximation 
of taking the gap to be (^-independent fails and large 
corrections are found. Indeed, at temperatures of order 
0.75T C , there is a twofold increase in the nonlinear co- 
efficient over that obtained with a ^-independent gap. 
This could have implications for the technological appli- 
cation of the cuprate materials as passive microwave de- 
vices for the communication industry. We also find that 
these corrections remain large in the presence of impurity 
scattering. 

A re- examination of the comparison of theory and data 
from the intermodulation power with this approximation 
removed finds that strong coupling effects would still 
be required to obtain a fit over the entire temperature 
range. However, this requires calculations of such sophis- 
tication that they are beyond the scope of this paper. 
The signature in support of d-wave symmetry remains 
clear, regardless. More importantly, a recent analysis of 
data from several films assumes that the explanation of 
the variation in the nonlinear coefficient and penetration 
depth results from impurity scattering. As the data is 
taken at high temperature, where the decay of the gap 
by the current is significant, we find upon reevaluation of 
the theory that the relationship between these two quan- 
tities can vary markedly and a different slope is obtained 
along with a different overall scale factor. Unfortunately, 
we must point out that, while a possible fit to the data 
still remains, the basic assumption that impurity scatter- 
ing is the source of the changes seen in the data implies a 
large variation in the T c of the different samples, which is 
not seen experimentally. Further experiments along these 
lines might help to resolve this issue more completely. 
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Finally, due to a recent experimental observation of 
possible d+s-wave gap symmetry in YBCO measured by 
angled-resolved electron tunneling, we have examined the 
nonlinear coefficient at both low T and T — > T c to de- 
termine signatures of the s-wave component. Indeed, at 
T — > T c , the isotropic behavior found for pure s- and 
pure d-wave is lost and a large anisotropy develops for 
d+s-wave as a function of the direction of the current 
in the plane. At low temperatures, the 1/T dependence 
remains but the fourfold anistropy of the coefficient is al- 
tered to twofold and the magnitude is modified. There- 
fore, should it become possible to do experiments where 



the direction of the current can be varied with respect to 
the antinodes or nodes, then we predict that signatures 
of a d+s-wave order parameter would be observable. 
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